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ABSTRACT 

Context. Chromospheres and coronae are common phenomena on solar-type stars. Understanding the energy transfer to these heated atmospheric 
layers requires direct access to the relevant empirical data. Study of these structures has, by and large, been limited to the Sun thus far 
Aims. The region of the temperature reversal can be directly observed only in the far infrared and submillimetre spectral regime. We aim at the 
determination of the characteristics of the atmosphere in the region of the temperature minimum of the solar sister star a Cen A. This will as a 
bonus also provide a detailed mapping of the spectral energy distribution, i.e. knowledge that is crucial when searching for faint, Kuiper belt-like 
dust emission around other stars. 

Metliods. For the nearby binary system a Cen, stellar parameters are known with high accuracy from measurements. For the basic model param- 
eters Teff. logg and [Fe/H], we interpolate in the grid of GAIA/PHOENIX stellar model atmospheres and compute the corresponding model for 
the G2 V star a Cen A. Comparison with photometric measurements shows excellent agreement between observed photospheric data in the optical 
and infrared. For longer wavelengths, the modelled spectral energy distribution is compared to 5pi/zer-MIPS, Herscliel-P ACS, Herscliel-SPIRE 
and APEX-LABOCA photometry. A specifically tailored Uppsala model based on the MARCS code and extending further in wavelength is used 
to gauge the emission characteristics of cfCen A in the far infared. 

Results. Similar to the Sun, the far infrared (FIR) emission of a Cen A originates in the minimum temperature region above the stellar photo- 
sphere in the visible. However, in comparison with the solar case, the FIR photosphere of a Cen A appears marginally cooler, r,^i„ ~ Ti(,o^^ = 
3920 ± 375 K. Beyond the minimum near 160 yum, the brightness temperatures increase and this radiation likely originates in warmer regions of 
the chromosphere of a Cen A. 

Conclusions. To the best of our knowledge this is the first time a temperature minimum has been directly measured on a main-sequence star other 
than the Sun. 
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1. Introduction 

aCentauri is the nearest stellar system to the Sun, located at a 
distance of 1.3 pc (tt = 747.1 ± 1 .2 mas, ,S6derhjelm, 1999 ). The 
physical binary is composed of two solar-like stars, the brighter 
of which, a Cen A (HIP 7 1 683, HP 128620) a G2 V star, is often 
considered a "solar twin" (jCayrel de Strobel|1996[|Melend ez eT| 
[ar|[2009). The companion, o-CenB, is of slightly later spectral 
type (Kl) and was recently found to host an Earth-mass planet 
( [Dumusque et al. 2012[ l. A possible link between chemical com- 
position in the atmospheres of solar twins and the formation of 
systems containing rocky planets was proposed by Melendez et 
|ar]( |2009l ). 

No planet has yet been found around the primary, however, 
but like the Sun, o-CenA shows evidence for chromospheric 



emission in the optical and ultraviolet spectral regions (Ayres 
et al.|1976 [Judge et aL 2004 and references therein) and should 
therefore also have atmospheiic regions where the temperature 
gradient turns from being negative to becoming positive. There, 
the radial temperature profile of the star should exhibit a mini- 



* Based on observations with Herschel which is an ESA space ob- 
servatory with science instruments provided by European-led Principal 
Investigator consortia and with important participation from NASA. 



mum. The rise in temperature beyond the "minimum tempera- 
ture" region is caused by non-radiative energy being deposited 
which leads to the heating of the higher atmospheric levels. The 
responsible physical processes are as yet un-identified and con- 
stitute the subject of intense study of solar physics and stellar 



astrophysics ([K alkofen 2007[ Wedemeyer-Bohm S. 
Cohen et al.|2005,,Haiper et al.|2012| i. 



etal. 2012 



The minimum temperature of the solar atmosphere can be 



measured directly only in the far infared ( Eddy et al.|1969[ Avrett 
2003 1. In the wavelength region ~ 50 - 350 fim, the atmosphere 
is becoming optically thick due to the dominating H~ free-free 



opacity (ocA , Geltman 1965 Doughty & Eraser 1966 1 and, con 



sequently, radiates at lower temperatures than the layers beneath 
(the photosphere in the visible, where To.s^m > !)■ The precise 
location of the temperature minimum is dependent on the de- 
tailed structure of the atmosphere and can, with no convincing 
theory at hand, be determined only from direct measurement. It 
is in this region, where non-radiative energy is deposited and its 
physics is of great general interest. 

By analogy, this phenomenon could also be expected to be 
found on a Cen A. We therefore set out to obtain such data with 
the space observatory Herschel (Pilbratt et al. |2010 l in the far 
infrared (FIR) and complementary ones with the ground-based 
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Table 1. Photometry and radiation temperatures of a Centauri A 



(jt/m) 


Sv 

(Jv) 




(K) 


Photometry 
& Reference 


0.440 


2215 : 


t41 


5792 


± 19 


B(l) 


0.550 


3640: 


t67 


5830 


±24 


V(l) 


0.790 


4814: 


t89 


5775 


±32 


1(1) 


0.440 


2356: 


1:43 


5856 


± 19 


B(2) 


0.550 


3606: 


t 66 


5818 


±23 


V(2) 


0.640 


4259: 


1:78 


5787 


±27 


R,,(2) 


0.790 


4784: 


1:88 


5764 


±32 


Ic (2) 


1.215 


4658: 


t 86 


5928 


±47 


J (3) 


1.654 


3744: 


t 69 


6121 


±60 


H(3) 


2.179 


2561 : 


t47 


5994 


±67 


K(3) 


3.547 


1194: 


t22 


5856 


± 78 


L(3) 


4.769 


592 ± 


11 


4995 


±69 


M(3) 


24 


28.53 


±0.58 


4736 


±91 


MIPS (4) 


70 
/U 




u. /u 


4599: 


t937 


Nllro ) 


70 


3.35 ± 


0.028 


4540 


±37 


PACS (5) 


100 


1.41 ± 


0.05 


3909: 


t 135 


PACS (6) 


160 


0.56 ± 


0.06 


3920: 


t394 


• PACS (5), (6) 


250 


0.24 ± 


0.05 


4084: 


t845 


•SPIRE (5) 


350 


0.145 


± 0.028 


4822: 


t927 


•SPIRE (5) 


500 


0.08 ± 


0.03 


5421 ± 


2018 


•SPIRE (5) 


870 


0.028 


± 0.007 


5738 ± 


1432 


•LABOCA (7) 



18" 



* Asterisks indicate values determined according to Eq. 1 . 
(1) HIPPARCOS, (2) Bes.sell ( 19901, (3) jEngels etal.| l [MT i. 

(4) G.Bryden [priv. com.; FWHM(24yum) = 6", (70//m) 
Binary separation on 9 April, 2005, 10'.'4. 

(5) Hi-GAL: KPOT_smolinar_l, fields 314_0 & 316_0. Herschel- 
beams FWHM(70yum) = 5'.'6, (100 /jm) = 6'.'8, (160 ^m) = 11'.'3, 
(250 yum) = 17'.'6, (350 yum) = 23'.'9, (500 yum) = 35'.'2. Binary sepa- 
ration on 21 August, 2010, 6'.'3. 

(6) DUNES: KPOT_ceiroa_l. Binary separation 29 July, 201 1, 5'.'7. 

(7) 384.C-1025, 380.C-3044(A): FWHM(870^m) = 19'.'5. Binary sep- 
aration 20-13 November, 2007, 8'.'8 and 19 September, 2009, 7'.'0. 



Table 2. Stellar parameters for a Cen A 



Ref^l)|1 
0(|2C 



Thevenin et al.'("2002 
2010jl. \Bonanno et al.,^2002j give Tq 



(2) Kervella et al.|p003l >, (3) Torres ef 
" - 4.57 ± 0.11 Gyr. 



APEX sub-millimetre (submm) telescope. These facilities al- 
low photometric imaging observations with high sensitivity. We 
present our findings in this Letter, which is organised as follows: 
in Sect. 2, the observations and reduction of the data are briefly 
described. The results are presented and discussed in Sect. 3, and 
in Sect. 4, our main conclusions are summarised. 



2. Observations and data reduction 



On July 29, 2011, PACS scan maps ( [Poglitsch et al.||20T0| of 
a Cen were obtained for the DUNES programme ( |Eiroa et al.| 



|2013[ ) at lOO/im and 160/im at two array orientations (70° and 
110°) to suppress detector striping. The selected scan speed 
was intermediate, i.e. 20"s"', determining the FWHM at the 
two wavelengths (6'.'8 and 11'.'3, respectively). In addition, 
PACS 70jum, 160jum and SPIRE 250/zm, 350 jum, 500 jum 
( [Griffin et aL]|2010|l data obtained as part of the Hi-GAL 
programme ( jMolinari et aL] |2010[ ) were also analysed (see 
http : //herschel . esac . esa . int/Docs/SPIRE/html/spire_oin . 
html#xl-980805.2.7). Our reduction work is thoroughly 



described by |Eiroa et al. (2013 i and Wiegert et al. (2013 1 
and the PACS calibration scheme is reported in detail in 
http : //herschel . esac . esa . int/twiki/bin/view/Public/Pacs 
CalibrationWeb#PACS_instrument_and_calibration. Aperture 
corrections of the flux and sky noise corrections for correlated 
noise in the super-sampled images were done according to the 
technical notes PICC-ME-TN-037 and PICC-ME-TN-038 and 
https : / /nhscsci . ipac . caltech . edu/ sc/index . php/Pacs/Ab- 
soluteCalibration. Both the relative and absolute positions in 
the sky are well known ( Pourbaix et al.|[2002| Table[l]l and are 
shown, for the epochs of the observations, in Fig. 1 of [Wiegert 
|etar] ( |20T3] l. 

The LABOCA ( Siringo et al. 2009 and http : //www . apex- 
telescope . org/telescope/) observations of a Cen were 
made during two runs, viz. in 2007, November 10-13, and in 
2009, September 19. The data associated with the programmes 
380.C-3044(A) and 384.C-1025(A) were retrieved from the 
ESO archive. The map data were reduced and calibrated using 
the software package CRUSH 2 developed by Attila Kovacs, see 
http : //www . submm . caltech . edu/~sharc/crush/download . htm. 
The data have been smoothed with a Gaussian of HPBW=13", 
resulting in an effective FWHM=23'.'4. However, fluxes in 
Jy/beam are given for an FWHM=19'.'5. 

All details regarding the data for both a Cen A and a Cen B 
will be presented in a forthcoming paper ( Wiegert et al.|2013"] l. 



Parameter, P (unit) 


Value of P 


i^P/P) 


Ref. 




Sun 


Q-Cen A 


Q-Cen A 




Mass, M (Mo) 


1.000 


1.105 


0.006 


(1) 


Radius, R {Rq) 


1.000 


1.224 


0.002 


(2) 


Radiative luminosity, L (Lq) 


1.000 


1.549 


0.042 


(3) 


Efi^ective temperature, T^ff (K) 


5770 


5824 


0.004 


(3) 


Surface gravity, log g (cgs) 


4.440 


4.306 


0.001 


(3) 


Metallicity, [Fe/H] 


0.000 


+0.24 


0.166 


(3) 


Age, r (Gyr) 


4.57+ 


4.85 


0.103 


(1) 



3. Results and discussion 

3. 1 . The spectral energy distribution of a Cen A 

Centauri was clearly detected at all observed wavelengths, 
with the possible exception at 500 yum because of high back- 
ground emission. Flux densities for a Cen A with their statisti- 
cal error estimates, from 0.4 to 870 fj.m, are provided in Table[T] 
together with the references to the photometry. For data where 
the binary was spatially unresolved (indicated by asterisks in the 
table), we used the PHOENIX models for both A and B to es- 



timate the relative flux contributed by a Cen A (Wiegert et al. 
|20T3l ), viz. 



> v,A 



> V, A 



+ S 



> 0.69 & < 1.0 if A < lOjum 
^ g I = const. = 0.69 if A> lOjum 



(1) 



This is in fair agreement with the PACS measurement at 70 yum, 
where the pair is resolved and the flux ratio Sa/Sb = 2.45 + 0.45, 
as compared to the PHOENIX value of 2.25. We are aware that 
using a constant value may not yield the correct answer at all 
wavelengths. However, with its deeper and more compact con- 
vection zone a Cen B is the more active and more X-ray lumi- 
nous star of the binary system ( DeWarf et al.|20l"0) l. Assuming 
that the wavelength of the temperature minimum is similar for 
both stars yields for a Cen A an optically determined Tj^m which 
is higher by 20% ( Ayres et al.|1976[ l. We will keep this caveat in 
mind when searching for the temperature minimum in a Cen A. 
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The observations of gCen are part of the DUNES pro- 3.2. The temperature minimum of a Cen A 



gramme (Eiroa et al. 2013 i which focusses on nearby solar 
type stars. The observations with Herscliel-PACS (Poglitsch et 
[al. 2010 ) at lOOyum and 160yum are aiming at detecting the stel- 
lar photospheres atanS /N > 5 at 100 ;um. Prior to Herschel, and 
surprisingly perhaps, not many data at long wavelengths and of 
high photometric quality are available for these stars. One rea- 
son likely being detector saturation issues due to their brightness 
(e.g., WISE bands W1-W4), another due to contamination in the 
large beams by confusing emission near the galactic plane (e.g., 
IRAS, ISO-PHOT and AKARI data). Another issue being how 
to define the photosphere at FIR wavelengths - which, in fact, is 
at the heart of this paper. 

Being as close-by as 1.3 pc and at comparable age, cCen A 
is an excellent astrophysical laboratory for stars very similar to 
the Sun. From numerous literature sources, [Torres et al.| ( |2010| l 
compiled the currently best available basic stellar parameters 
and the estimated errors on the physical quantities are gener- 
ally small (Table|2|. However, whereas the tabulated uncertainty 
of the effective temperature of a Cen A is less than half a per- 
cent, the observed spread in Table 1 of |Porto de Mello et al. 
( 2008| l does correspond to more than ten times this much. On 
the other hand, the radius given by Torres et ajj ( |2010| l is that 
directly measured by Kervella et al. ( 2003 ) using interferometry 
(and corrected for limb darkening), with an error of 0.2% ( |Bigot| 
|et al.|20()6| . The mass has been obtained from asteroseismolog- 
ical measurements and is good to within 0.6% (Thevenin et al. 
2002[ l. However, it is also evident from the table that the twinship 



of a Cen A and the Sun is not identically close in every detail. 

For such an impressive record of accuracy for the stellar pa- 
rameters of the a Cen components it should be possible to con- 
struct theoretical model photospheres with which observations 
can be directly compared to a high level of precision. In this pa- 
per, we report on FIR and submm observations, which should 
provide valuable constraints on the Spectral Energy Distribution 
(SED). This could potentially be useful to gauge the tempera- 
ture minimum at the base of the stellar chromosphere. A clear 
understanding of the latter is crucial when attempting to deter- 
mine extremely low levels of cool circumstellar dust emission 
( [Eiroa et aI.||20II[ |20I3) l. Because of binary dynamics, dust is 
not expected to contribute to the stellar radiation from a Cen A 
( [Wiegert et al.|2013, and references therein). 

The atmosphere model for a Cen A has been computed by 
a 3D interpolation in the high-resolution PHOENIX grid for 
GAIA (Brott & Hauschildt ,2005 ) and with the following pa- 
rameters: (reff, \ogg, [Fe/H])=(5824K, 4.306, +0.24), where 
[Fe/H]* = log(A?Fe/A^H)* - log(A?Fe/A^H)Q. This model is shown 
in Fig.fT] together with the excellently agreeing photometry 
(TablefTT The photometry has not been used in the analysis, but 
serves to illustrate the good agreement between this model and 
the observations. 

The PHOENIX model spectra are computed up to 45 ;um. 
These models are used extensively by the DUNES programme 
and therefore also exploited here. However, for this particu- 
lar study, we have also used a specifically tailored model (the 
"Uppsala model") based on the MARCS code (Gustafsson et al. 
[2008| l, with, of course, exactly the same atmosphere parameters 
as before. This model computation extends to atmospheric lay- 
ers more than 1500 km up, including those regions which emit 
predominantly around 200 jum. The infrared part of the Uppsala 
model is shown superposed onto the PHOENIX model in Fig.[T] 
where it can be seen that these two atmosphere models are vir- 
tually indistinguishable. 



The observed FIR fluxes of a Cen A seem somewhat lower than 
in the model but appear to turn upward in the submm/mm bands 
(inset of Fig.[T). The ID Uppsala (or PHOENIX) model de- 
scribes stellar atmospheres in local thermodynamic equilibrium 
(LTE) and does not account for any temperature inversions. Both 
non-LTE effects and the change of the temperature gradient will 



influence the emergent spectrum (for details, see De la Luz et 
al. 2011). In the far infrared, this emission from the solar chro- 
mosphere exhibits a minimum in brightness temperature around 



150jum (Avrett 



2003 and references therein). A chromosphere 



has also been observed for a Cen A in ultraviolet line emission 
(e.g.. Judge et al.|2004[ l. From the interpretation of the wings of 



the optical Call K line, Ayres et aL] (1976 ) deduced a value of 
TmmlTeS - 0.78, similar to the solar value of 0.77. According 
to Judge et al. ( 2004]), the level of activity of a Cen A is low and 
similar to that of the Sun in an intermediate state of its cycle and, 
on the basis of long-time X-ray monitoring, this apparent lack of 
variability was also emphasised by Ay res (2009 ). 

We assume an atmosphere in radiative and hydrostatic equi- 
librium. The flux received at the Earth, Sy, is related to the out- 
ward flux through the surface of the star, Fy (e.g., ,Mihalas 1978) , 
through 



'(Iff 



(2) 



where (py is the angular diameter of the stellar disc of radius 
Ry, at a given frequency v or wavelength A, and as seen at a dis- 
tance D. For these parallel stellar rays, a brightness temperature, 
rB(v), can be defined through the Planck function by 



Fy^nByin) , 



so that with 0v/2 = RylD, 



Tb(v) 



h V 
~k 



(3) 



(4) 



In the far infrared, /?fir = ^0.5 + AT?, where /?o 5 is the "photo- 
spheric" radius of the star. There the optical depth in the visible 
is unity (tq.s ~ 1). The temperature minimum is found in regions 
higher up, where tfir > 1 (but T0.5 <K 1). It is straightforward to 
show (e.g., Tatebe & Townes|2006[ l that AR/Rpi^ is of the order 
of a few times 10 Therefore, AR corresponds to some 500 km, 
and will per se not introduce any significant errors for stars like 
the Sun and a Cen A (luminosity class V) and in Eq. 3, we use 
^FiR - R0.5 (cf. Table|2]i. However, over this distance, the density 
drops by an order of magnitude, making model computations in- 
creasingly difficult. 

In Fig.[T| T^iv) is shown for the Uppsala model photosphere, 
where LTE is assumed for the free-free H" continuum, together 
with a semi-empirical chromospheric model of the quiet Sun 
(VALIIIC, Vernazza et al.|1 981 ). Solar data are shown as open 
circles and the observations of a Cen A as filled squares. 

Using the Uppsala model, we find that;i'ioo/jm = -4.4, where 
Xy - (5v,obs - 'Sv,mod)/o"v. The difference in brightness tempera- 
ture by 500 K in the 100 to 160 fim region could therefore be sig- 
nificant and it cannot be excluded that in a Cen A the atmosphere 
becomes optically thick at higher, cooler levels than the model 
and that the structure is different from that in the SurQ Important 

' Actually, a very similar scenario is proposed also for the Sun by 
|Ayres[ j2002[ l in order to explain the lower T,^],, indicated by observa- 
tions of infrared CO lines. 
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is however the fact that the FIR-SED indeed goes through a min- 
imum, as the brightness temperatures rise at longer wavelengths 
and, at e.g. 870//m, H free-free emission from the stellar chro- 
mosphere appears to dominate. 

In the solar atmosphere, the temperature minimum occurs 
around 150 jum. This seems to be the case also for aCenA, 
where ri(,o/jm = 3920 ±375 K. It is customary to express the min- 
imum temperature of the stellar atmosphere in terms of the effec- 
tive tempera t ure, Te ff, and for aCenA, T^in/Teff = 0.67 + 0.06. 
Ayres et al. ( 19761) had earlier, from Call K-line fitting, esti- 
mated this ratio as 0.78 to 0.79, for an assumed T^ff - 5770 
to 5700 K, hence T^i^ ~ 4500 K, which would indicate an "op- 
tical minimum temperature" about 500 K higher than what we 
determined from a direct measurement in the far infrared. This 
result may seem surprising, but we have to recall that similar is 
seen in the Sun, where diff'erent diagnostics yield a span in min- 
imum temperature of about 600 K ( Avrett 2003). It is also worth 
remembering that even for the Sun, a unique model for its chro- 
mospheric emission has yet to be found (e.g., |De la Luz et al. 
|201 I V For its "sister star" ffCenA, future data at submm and 
mm wavelengths, with for instance ALMA, should contribute to 
a better characterisation of the atmosphere of our nearest neigh- 
bour. 



4. Conclusions 

• We successfully observed the far infrared energy distribu- 
tion of the nearby G2 V star a Cen A with instruments aboard 
Herschel and from the ground. 

• The observed radiation temperatures appear lower than what 
is expected on the basis of extensions of LTE stellar at- 
mosphere models. Near 160 yum, the minimum temperature 
in the atmosphere of a Cen A, T„^i„/Tf.ff - 0.67 + 0.06, is 
marginally lower than the ratio of 0.73 observed in the Sun 
at 150;um. 

• At 870 /vm, however, the emission from a Cen A appears to 
originate from regions at higher temperatures, as might be 
expected, viz. T^jo^i^/Tstf = 1.0 + 0.2. 

• Cold temperature minimum regions will lead to the underes- 
timation of the amount of dust present in cold debris discs. 
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Fig. 1. The SED of a Cen A, where the blue line represents the PHOENIX model photosphere (computed up to /I = 45 /im and 
Rayleigh- Jeans extrapolated beyond) and the thick (purple) line the Uppsala model photosphere (extending to A - 200 yum). 
Photometric data are shown with their Icr error bar estimates (Table[T]l. The inset shows the brightness temperature of a Cen A 
(squares) and the Quiet Sun (circles) in the far infrared and submm. The solar data from the compilations by jGuetal. ( 1997| ) and 
Loukitcheva et al.| ( )200 4) are shown together with a semi-empirical chromosphere model for the Sun ( Vernazza et al.|198T black 
dashes: VALIIIC). At 70jum, the symbol with the larger error bar (green) represents the Sp/fzer-MlPS datum and filled (blue) 
squares i/eric/zeZ-PACS/SPIRE and LABOCA data. Horisontal bars indicate filter widths. 



5 



